A Pinned Polymer Model of Posture Control 

Carson C. Chow 1 and James J. Collins 1 ' 2 
1 NeuroMuscular Research Center and 2 Department of Biomedical Engineering 
Boston University, 44 Cummington Street, Boston, Massachusetts 02215 

(January 11, 1995) 

Abstract 

A phenomenological model of human posture control is posited. The dy- 
namics are modelled as an elastically pinned polymer under the influence of 
noise. The model accurately reproduces the two-point correlation function of 
experimental posture data and makes predictions for the response function of 
the postural control system. The physiological and clinical significance of the 
model is discussed. 

PACS numbers: 87.45. Dr, 0.5.40.+J 
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I. INTRODUCTION 



The human postural control system is highly evolved and complex. Recent analyses of 
quiet-standing posture data suggest that stochastically driven dynamics may be present in 
the system [1,2]. Here we present a phenomenological Langevin equation that captures the 
underlying physical mechanisms and reproduces the experimental posture data. The model 
is simple enough to allow quantitative analyses. Here we obtain correlation and response 
functions for the model and relate the properties of these functions back to physiological 
origins. 

In Ref. [1], Collins and De Luca measured and analyzed the time- varying displacements 
of the center of pressure (COP) under the feet of quietly standing subjects. They found 
that the dynamics of the COP time series were not consistent with low-dimensional chaos. 
Subsequently, they constructed the two-point correlation function in the anteroposterior 
(front-to-back) direction, j/, defined by C(t 2 — ti) =< ( y (t 2 ) — y(ti)) 2 >, from the data. An 
example is shown in Fig. 1. Collins and De Luca found that for healthy young subjects 
the COP two-point correlation function exhibited three approximate scaling regions where 
C(t) ~ t 2H [3]. In the first region H 4/5, in the second H 1/4, and in the third 
H 0. These regions can be described as: (1) drift-like or free streaming-like, (2) diffusive, 
and (3) bounded or saturated [4]. Together the above posture results suggest a continuum 
phenomenological stochastic model for posture control. 



We consider the transverse motion of the human body in one spatial dimension, y. We 
would like to model the dynamics of this motion as a function of the height variable, z, and 
time, t. The simplest continuum model of the body is that of a flexible string or polymer 
under tension with friction 



II. THE MODEL 



pd 2 t y + fid t y = Td 2 y, 
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where p is the mass density, p is a friction coefficient and T is the tension. In upright 
stance, the body has an inverted pendulum-like instability which is countered by a control 
system. The precise form of this control system is highly complicated and depends on 
many physiological processes [5]. A simple approximation is to consider the body to be 
"pinned" to the upright position by an elastic force. In terms of the model, this would 
amount to embedding the polymer in an elastic sheet. The corrections of the control system 
cannot always be perfect [6] and this is modelled by a stochastic force. The entire postural 
control system is reduced to an elastically pinned polymer under the influence of stochastic 
fluctuations. The full equation takes the form 



where K is the elastic restoring constant of the pinning and F is the stochastic driving force. 

Equation (2.2) can be put into a more convenient form for analysis by dividing through 
by fi and relabelling the constants: 



In this form, f3 and a -1 have dimensions of time, and v has the dimensions of length squared 
divided by time. By rescaling the length, the parameter v can also be eliminated; however, 
in order to keep length and time dimensions explicit, it is retained. The free parameters 
in the system are the two time scales f3 and a -1 and those associated with the stochastic 
forcing. 

The stochastic forcing rj(z,t) is assumed to have a correlator of 



It is uncorrelated along z but can possess temporal correlations represented by the function 
D(t). In Fourier space we assume the correlator can be approximated by the form 
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This assumes that for long times (t » <5 _1 ) the forcing is uncorrelated but for short times it 
has some power-law scaling, < ^(lo 1 )^^) >oc u~ 2l 8{u}' + u). These short-time correlations 
are assumed to arise from neuromuscular effects, such as intrinsic muscle force fluctuations 
and/or feedback-loop delays. The values of the parameters will be constrained by the data. 

The dynamics of the model are to be compared to the measured motion of the COP. 
It is assumed that this motion can be effectively mimicked by the motion of a particular 
point along the polymer (body). Hence, it is necessary to compute the two-point correlation 
function C(t 2 — ti) =< (y(z 0} t 2 ) — y(z 0} ti)) 2 > taken at a specific point z = z . In actual 
posture control for quiet standing, the feet are essentially fixed to the floor and hence the 
pinning strength is infinite there. As one goes upwards, the pinning strength is reduced. 
The analysis of a model with inhomogeneous pinning, however, is difficult. Nonetheless, 
since only the motion at one point along the body is required, we can consider a model 
with an "averaged" constant pinning strength. We also assume that the boundary effects 
are unimportant and consider an infinitely long polymer. In this case, every point along the 
polymer will have equivalent dynamics. These approximations are valid for strong pinning 
and for the dynamics near the center of mass away from the boundaries which we assume is 
a location that effectively mimics the COP dynamics. 

The correlation function can be computed from Eq. (2.3). It is most convenient to 
analyze the problem in the Fourier domain. Fourier transforming Eq. (2.3) in space and 
time and solving for y(k } io) yields 





The auto-correlation function is then formed to yield 





where Eq. (2.5) has been applied. 



Inverse Fourier transforming in W 



'', k' and k yields 



<y(z,u)y(z',u) >= 2D6 2 ^ [ 



dk 




'f 2tt [(/3uj 2 - vk 2 - a) 2 + lo 2 ]{lo 2 + S 2 )^ ' 
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where the integration is along the contour F . We want the single-point temporal correlations 
so we set z = z' = z in Eq. (2.8) and obtain 

DS 21 f 1 

S(u) =<y(z ,u))y(z ,u))>= / — — ^- c2\~ t dk - ( 2 - 9 ) 

7r jf ((pur — z//r — a)^ + uj l ){uj l + o^) 7 

The two-point correlation function is given by 

C(r) = 2(5(0) -5(r)) (2.10) 

where S'(t) is the temporal auto-correlation function obtained by inverse Fourier transform- 
ing Eq. (2.9). 

The different scaling regions in frequency and hence time can be obtained by taking 
suitable limits of Eq. (2.9). Note that vk 2 /u>, f3u and a/u are dimensionless quantities. 
With this in mind Eq. (2.9) can be rewritten as 

D 1 r S 21 

S ^ = ~Vuj 2 (uj 2 + 8 2 ^ If ((J3uj-a/uj -uk 2 /uj) 2 + l) dk ' 

Consider the range of parameters where p 1-1 ~ 8 » u » a. Then Eq. (2.11) can be 

approximated by 

~ if 3/ 2 r -TTT du > ( 2 - 12 ) 



where u = ^Jv/tok. The integral is convergent so in this regime S(to) oc lo~ 3 / 2 . By dimen- 
sional analysis, this translates to a scaling behavior of S(t) ~ r 1 / 2 , in the time domain [7]. 
This corresponds to the middle scaling region (where H 1/4) in the COP data of Fig. 1. 
This scaling region is valid in the range f3 << r << a -1 . 

Now consider the long-time region where u « a « p 1-1 ~ 8. Equation (2.9) then 
takes the form 

S(u) ~ - r 1 dk. (2.13) 
7r J-oo \vk z + ay 

The integral in Eq. (2.13) is convergent so we find that S(u) oc u>°, which gives S(t) ~ 
for r >> a -1 . This corresponds to the saturated region (where H 0) in the COP data in 
Fig. 1. 



Finally, in the short-time region u » f3 1 >> 1, <5 ~ /3 we have u 2 >> a and u 2 f3 >> 1 
which leads to 

DS 21 f°° 1 

s{LO) ~ ^/^/^ (1 _„ 2)2 + £ > ( 2 - 14 ) 

where 



(2.15) 



/3 ftto 

In the limit — > oo (e — > 0), the integral in Eq. (2.14) is proportional to e -1 so S(u) behaves 
as u~ 2 ( 1+1 \ By dimensional arguments we find in the time domain, Sir) ~ t^ 1+1 \ 

The behavior of the complete S{u) is shown in Fig. 2. The calculated two-point correla- 
tion function takes the form corresponding to the COP data in Fig. 1. Fitting the parameters 
to the data in Fig. 1 yields values of a -1 10 s, f3 ~ 0.5 s, and 7 0.72. 

The saturation width C(r = 00) can also be estimated. This value represents the 
maximum transverse extent the body attains during quiet standing. From Eq. (2.10) we see 
that C(oo) = 2S(t = 0) since S(t = 00) = 0. Hence C(oo) = 2S(t = 0) = 2 Jf^ S{to)dto. 
For simplicity consider 7 = 0, i.e., no correlations in the noise. The temporal correlations 
in the system noise in healthy subjects only extend over short time periods and do not play 
a role in the saturation width. By making the substitution u = \J f3 / au and v = \jv j ' ak in 
Eq. (2.9) we obtain 

0) ^ f°° f°° ^ U ^' V (2 16) 

a\fv~ft J-00 J-00 (u 2 — v 2 — l) 2 + u 2 /af3 

D f°° dv , 

(2.17) 



[V 2 + l) 3 / 2 ' 

i.e., C(oo) oc D / y/ pa. By absorbing the constants into D } an effective noise amplitude 
D e {t = C(oo)\/a can be defined. (Recall that the tension parameter v is not independent 
and could have been scaled out into the noise amplitude from the outset.) From the COP 
data in Fig. 1, we get D e ff 60 mm 2 s~ 1/ ' 2 . 

The values of the respective parameters have implications for the postural control system. 
The non-zero value of 7 implies short-time correlations in the noise. This corroborates the 
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open-loop control hypothesis in Refs [1,2]. In these short times, the "kicks" produced by 
the system are correlated to each other and are not reacting to the position of the COP. 
The body responds inertially to these kicks, resulting in drift-like dynamics. The second 
scaling region of H 1/4 represents "diffusive" scaling. In this region, the dynamics involve 
a balance between the damping, tension, and uncorrelated noise of the system. Physically 
this implies that the body is aligning itself independently of whether or not it is upright. In 
the saturated region, the pinning is finally felt and the body then attempts to return itself 
to an upright position. The noise amplitude does not affect the scaling regions, but it does 
establish the saturation width. The physiological significance and relevance of these and 
different parameters are discussed in Sec. IV. 



III. RESPONSE TO PERTURBATIONS 

With the pinned polymer model, the response to a perturbing force can also be calcu- 
lated. From a physiological standpoint, this is important given that a significant number 
of contemporary investigations in posture control have involved analyses of the measured 
response of the human body to various external perturbations [8]. This issue is discussed 
further in Sec. IV. 

To obtain the response function consider the polymer model Eq. (2.3) with an imposed 
perturbation 

f3d 2 y + d t y = vd 2 z y - ay + r)(z, t) + e(z, t), (3.1) 

where e(z,t) is a spatiotemporal perturbation. For simplicity we consider 

e(z,t) = e6(z)6(t). (3.2) 

This corresponds to an impulsive kick at z = 0. We Fourier-Laplace transform Eq. (3.1) 
and obtain 

—uj 2 f3y — iujy = —uk 2 y — ay + fj(u!, k) + e. (3-3) 
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The response function is defined as 



g(k,u) = r-f) = — — . (3.4) 

v ; \8t -uj 2 f3 - iuj + uk 2 + a v ; 



Inverse transforming leads to 



f dco ■ , f°° dk 1 
g(z,t)= / —e- twt / — e lkt — : . 3.5 

yv ' ; Jl 2tt J -co 2tt -lo 2 (3 - iuj + vk 2 + a v ; 



The k integral can be performed to yield 



1 /■ (-l3w 2 -iw+a)/i/ 

9(z,t) = —= / due-*** 2 , ■ (3.6) 

By making the transformation u — > u — i/2/3 } the response function takes the form 



t/2/3 ,. e -zV^ 



2_,„2 



g(z,t) = - ir - F = / due~^ f , (3.7) 



where 

a = ^^, Z = ^-|*|. (3.8) 

For 4a/3 > 1, there are branch points at u = ±a. By definition there is no response before 
t < Z so the Laplace contour must be above the branch cut from —a to a. For t > Z } we 
can close the contour in the lower half plane and collapse the contour around the branch 
cut. The path of integration is then taken along the two sides of the cut. The integral is 
completed with the substitution u = a sin <f> and some algebraic manipulations to yield 

-t/2/3 i f27T 

g( z ,t) = e(t-Z)—=— e ia ^^^d<f>, (3.9) 
l\Jvp Air Jo 

where 0(t) is the step function. The integral in Eq. (3.9) is the integral representation for 
the zeroth-order Bessel function so 

p -*/2/3 , 

g (z,t) = 9(t - Z)^=J (aVt 2 - Z 2 ). (3.10) 



20^ 

For the case where 4a/3 < 1, the branch points lie along the imaginary axis. A rotation in 
the complex plane by tt/2 in Eq. (3.7) leads to the result 
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p -*/2/3 , 

g(z,t) = 9(t - Z)—=I (\a\VW=z5), (3.11) 

L\JVp 

where I is the zeroth-order modified Bessel function. At the critical value of Aa/3 = 1, the 
response is given by 

e -t/2/3 

g (z,t) = e(t-Z)—=. (3.12) 

l\Jvp 

If the response is measured at the location of the perturbing kick (z = 0), then we can 
set z = in the above expressions. The response for different values of (3 and a for z = are 
shown in Fig. 3. For small values of a, the response time is on the order of 2/3. If (3 is very 
small, \a\ becomes very large and the asymptotic form of the modified Bessel function can be 
used, i.e., Iq(x) ~ e x /\ // 2ttx for x — > oo. This then leads to the expected damped-diffusive 
response behavior of g(0,t) ~ e~ at j\J Airpt. For near-zero pinning (a ~ 0), we obtain a 
power-law response behavior of g(0,t) ~ t -1 / 2 . 



IV. DISCUSSION 

With the pinned polymer model, posture control is characterized in terms of four pa- 
rameters: (1) the onset-of-damping time scale /3, (2) the onset-of-pinning time scale a -1 , 
(3) the short-time noise correlation exponent 7, and (4) the effective noise amplitude -D e ff- 
These parameters could be used to form the basis of a physiological and clinical classifica- 
tion scheme. For instance, deviations in a parameter or set of parameters from those for 
a population of age-matched, healthy individuals could serve to indicate a possible balance 
disorder. Each parameter has some physiological implication and so the model can make 
predictions of how various physiological changes (e.g., due to injury or disease) would alter 
the COP data. 

For the quiet-standing posture data of healthy subjects [1,2], nonlinearities are not nec- 
essary — it appears that the pinning term is strong enough to negate any nonlinear effects. 
However, this may not be the case for certain patients with balance disorders. For instance, 
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nonlinearities will certainly play a role for large transverse excursions. Obvious effects in- 
clude stepping or falling. These catastrophic effects lie outside the realm of the model. Weak 
nonlinear effects could be incorporated by allowing the pinning coefficient a to be a function 
of y. This would be manifested as terms of the form y n in Eq. (2.3). The y — > —y symmetry 
of Eq. (2.3) could be broken for even n. This would make sense from a biomechanical 
standpoint since anteroposterior movements are intrinsically asymmetric [9]. 

The quiet-standing posture data [1,2] showed that short-term correlated noise is present 
in the postural control system, i.e., for short times the kicks produced by the system are 
temporally correlated. This finding was confirmed by the model, i.e., the stochastic forcing 
needed to have some power-law scaling over short times in order for the model's output 
to reproduce the experimental posture data. These results imply that during undisturbed 
stance the postural feedback mechanisms are not necessarily operational over short periods 
of time (t < 1 s). This is in line with the open-loop control hypothesis of Collins and De 
Luca [1,2]. Physiologically, this scheme could result of inherent noise, feedback-loop 

delays, feedback thresholds, and/or sensory "dead-zones". 

A discrete version of this model would involve a chain of coupled and bounded random 
walkers with short-time correlations. The coupling between the walkers is required for 
the diffusive region. The pinned polymer model considers an infinite chain but a finite 
number of coupled random walkers would probably suffice. The correlations could arise 
from time-delays in the stepping probabilities of the random walkers [10]. The advantage of 
a continuum version is the ease in which calculations can be made. 

Various behaviors could arise for different parameter regimes. For instance, if the damp- 
ing and pinning time scales were roughly equal (/3 ~ « _1 ), then the diffusive region would 
disappear. The COP would essentially drift until it was corrected by the elastic pinning. 
It is possible that a control strategy involving such dynamics would be advantageous for 
individuals who require quick movements (e.g., athletes) since the dynamics imply that the 
damping of the overall system is relatively weak. An individual utilizing such a strategy 
would be metastable and capable of making quick adjustments. 
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From the standpoint of the model, balance problems could be manifested in several ways. 
For instance, if the pinning strength were small, the COP would fluctuate diffusively until 
its transverse displacement exceeded the stability (balance-maintenance) threshold — this 
could result in a fall. If the damping were also weak so that f3 > a -1 , then oscillations 
would occur before saturation. Falls could be triggered if the amplitude of these oscillations 
was large enough to exceed the stability threshold. Similarly, if the noise amplitude was 
sufficiently large, then the amplitude of saturation could exceed the stability threshold — as 
with the aforementioned scenarios, this effect could trigger a fall. From the above discussion, 
it is clear that the pinned polymer model may be useful for extracting clinically relevant 
information from quiet-standing COP time series. 

The pinned polymer model has significant utility for understanding posture control quan- 
titatively. As noted above, the model can completely characterize the statistical behavior 
of quiet-standing COP trajectories in terms of four parameters. Once these parameters are 
fixed from the "quasi-static" posture data, the model can be used to make a prediction 
of the dynamic response of the postural control system. The model thus enables a direct 
comparison between the "quasi-static" and "dynamic" postural control systems. From a 
physiological standpoint, this is important because it allows one to test the hypothesis that 
the two control modes are equivalent, i.e., under quiet-standing and dynamic conditions, 
the postural control system utilizes the same mechanisms. This hypothesis has not been 
addressed previously because of the lack of a suitable quantitative model or approach. 

Experiments are under way to measure the dynamic response of the human postural 
control system. If the parameters from the quiet-standing posture data match those of the 
response function, this would support the aforementioned hypothesis and further validate 
the pinned polymer model. It is expected that the dynamic response will delineate when 
the two mechanisms differ. For large perturbing kicks, it is possible that nonlinear effects 
will appear and cause a deviation from linear response. 

The model does not address the actual neuromuscular control mechanisms involved in 
balance regulation. For instance, it is unknown whether the stochastic forcing is intrinsic in 
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the output of skeletal muscles or whether it arises from a combination of incomplete sensory 
information, feedback-loop delays, and/or feedback thresholds. In order to address these 
issues, a detailed model of posture control that accounts for neurophysiological components 
and body mechanics is required. However, the pinned polymer model does allow a quanti- 
tative characterization of posture control and an assessment of physiological consequences. 
Moreover, any detailed model must reduce to the pinned polymer model at a coarse-grained 
level. 
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FIG. 1. A resultant log-log plot (solid line) of the COP two-point correlation function (C(r), 
where r is time interval) for a healthy young subject. This plot was generated from five 90-s COP 
time series. Shown also are the fitted regression lines (dashed lines) and the computed values of 
the scaling exponents (H) for the respective scaling regions. 
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FIG. 3. Three example response functions corresponding to the parameters: v = 0.25, (3 = 1 
and a = 0.5 (dotted line), a = 0.25 (solid line), a = 0.128 (dashed line). 
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